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Abstract. In a previous paper it has been shown that the interference of the first and second order pole 
of the Green’s function at an exceptional point, as well as the interference of the first order poles in the 
vicinity of the exceptional point, gives rise to asymmetric scattering cross section profiles. In the present 
paper we demonstrate that these line profiles are indeed well described by the Beutler-Fano formula, and 
thus are genuine Fano resonances. Also further away from the exceptional points excellent agreement can 
be found by introducing energy dependent Fano parameters. 
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1 Introduction 

Fano resonances were first discovered in atomic physics [l 
|2j. In cross sections in nuclear physics they also appear anc 
are called Feshbach resonances |3][4]. The characteristic 
features of these resonances are asymmetric line shapes, 
caused by the interaction between a discrete state and a 
continuum of states, as described by the formula derived 
by Fano in his landmark work in 1961 5|. Later Fano 
himself developed the subject further in numerous papers, 
review articles [6], and a book 17 . 

To date Fano resonances have been observed in a wide 
array of different subfields of physics. They have been 
found, to give just a few examples, in mesoscopic con¬ 
densed matter systems [§], in light scattering by finite ob- 

in wave transmis- 


previous paper 21 . The present paper implements this 


stacles [9], in nanoscale structures 10 
sion through resonant scattering structures 11 , in broad¬ 


band atto-second-pulsed XUV liglit_ experiments 12 13 
and in coupled plasmonic systems 


14-17 


19-21 and references quoted therein, see also, e. g., 22- 


claim with substantial numerical evidence. 

We demonstrate that both at and in the vicinity of 
the exceptional point the cross sections indeed give rise 
to Fano profiles locally, and therefore can be interpreted 
as genuine Fano resonances. Moreover, the cross sections 
can globally be described by the Beutler-Fano formula, 
if energy-dependent Fano parameters are introduced. As 
an aside we demonstrate that specific parametrisations 
traditionally used |25| are bound to fail at the EP, since 


the two peaks of the cross section cannot be associated 
with individual resonances. 


2 The model 

For the reader’s convenience we briefly review the model 
introduced in 21 . The model is inspired by work of Joe et 


al. 26 , who had pointed out an analogy between quantum 


In this paper we take a different theoretical perspec¬ 
tive on Fano resonances. Similar to the well-established 
description of an isolated resonance by a pole in the com¬ 
plex plane of the scattering function, we here propose an 
understanding of Fano resonances by specific singularities 
denoted as exceptional points (EPs) 18]. 

There is a host of literature about exceptional points 


interference in Fano resonances and classical resonances in 
coupled oscillators (see also 27 28 ). The system consists 


24 . In our context the important point is the occurrence 
of an additional pole of second order in the scattering 
function at the EP. As a consequence, the cross section 
(the modulus squared of the scattering function) is then 
characterised by the all important interference of the poles 
of first and second order. This has been hinted at in a 


of two one-dimensional harmonic oscillators with unper¬ 
turbed eigenfrequencies uq, to 2 and damping constants k \, 
k‘ 2 , coupled by a spring with spring constant / and damp¬ 
ing constant g , and periodically driven by an external force 
with frequency w. We note that in the literature the model 
is used to describe the plasmonic response of coupled plas¬ 
monic structures (cf. 16|[l7 ). 

After setting up the equations of motion in phase space, 
(p, q) T = M(p, q) T + cexp(iwf), resonant solutions with 
real frequency w are found by the roots of the characteris¬ 
tic polynomial of the matrix A4. Additional singularites, 
called exceptional points, occur for complex values of u 
where in addition the first derivative of the characteristic 
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polynomial vanishes. This leads to 4 real-valued equations 
for the 8 real parameters uq, aq, fei, k 2 , g, /, 5J(w), 3(w). 
If we provide specific values for the eigenfrequencies and 
damping constants, we are left with 4 equations for the 
4 unknowns g, which in general have to be 

solved numerically. However, in the special case that the 
individual oscillators are undamped, k\ = /c 2 = 0, as used 
in the following, the exceptional points can be given in 
analytical form: 


3?[wep] 


S[wep] 


1 /(3w 2 + a*!) (w 2 + 3w|) 

2-\/2 V + uj % 

2\[2 \Jw\ + u>2 


9 ep = -S[wep] = 




2a/2 


/ep = 


4 (wf + • 


In close vicinity of an exceptional point the higher di¬ 
mensional eigenvalue problem can be reduced to an effec¬ 
tive two-channel scattering problem. The advantage of the 
reduction is the analytic availability of the Green’s func¬ 
tion and thus the scattering amplitudes both near and at 
the exceptional point. /From the Green’s function the T 

matrix can be calculated, and from it the cross sections 
|2 


|Tn| , IX 22 I 2 . For details we refer the reader to 


21 


3 Fano profiles 


Table 1. Resonance energies Eb, widths E, and asymmetry 
parameter q of the Fano profiles around the exceptional point 
shown in Fig. [2] 



Eb 

r 

Q 

a) 

2.04628 

0.00078 

-6.43165 

b) 

2.04242 

0.00289 

-3.03946 

c) 

2.03899 

0.00623 

-2.09385 

d) 

2.04917 

0.02755 

-0.06180 




energy E 


energy E 




energy E 


energy E 


Fig. 1. Fano fits (red) of the cross sections at and near the 
exceptional point. The parameters are on = 2.00 and uj 2 = 
2.10, ki = &2 = 0, the exceptional point appears at /ep = 
0.005, g E p = 0.0499, and ca B p = 2.0500 - 0.04999L (a) / = 
/ep + 1.0, (b) / = /ep + 0.5, (c) / = /ep + 0.3, (d) / = /ep. 


3.1 Fits with one Fano profile 


To prove that the asymmetric line profiles found in 21 


around exceptional points are actual Fano resonances, we 
first recall the Beutler-Fano formula (see e.g. |7|) 


o-(e) 


(e + g ) 2 

e 2 + 1 


(1) 


It quite generally describes line profiles for a physical pro¬ 
cess where a continuum state interacts with a bound state 
embedded in the continuum. In ([l]) 


E — Eb 

r /2 


( 2 ) 


is the reduced energy which measures the energy relative 
to the position of the resonance Eb in units of the half¬ 
width r/2 being the width of the quasi-bound state in the 
continuum (note that the interaction endows the quasi¬ 
bound state with a width). The parameter q determines 
the shape of the resonance and depends on the ratio of 
the transition matrix elements linking the initial state to 
the discrete and continuum parts of the final state. 

In the following we assume no damping of the individ¬ 
ual oscillators and choose uq = 2.00 and uq = 2.10. The 


exceptional point then appears at /ep = 0.005, gEP = 
0.0499, and wep = 2.0500 — 0.04999L In Fig. |T]we illus¬ 
trate Fano fits of the cross sections at and near the excep¬ 
tional point. The fitting parameters are Eb , F, and q , in 
addition we have a scaling factor in front of eq. <[T|) to scale 
the height of the resonance and a global shift of the cross 
section to account for a background. Numerical values of 
the fitting parameters are given in Table [l] 

It is evident from Fig. [l] that both at and further away 
from the exceptional point the Beutler-Fano formula lo¬ 
cally describe the asymmetric line profiles perfectly. It 
thus confirms that they are Fano resonances. 


3.2 Fits with two Fano parameters 

If one allows for two sets of energy dependent Fano pa¬ 
rameters the two peaks in the cross sections and the asym¬ 
metric line profile can indeed be modelled over a far wider 
range of the energy. Examples are shown in Fig. [2j the 
fitting parameters are listed in Table [2] It can be seen 
that further away from the exceptional point the profiles 
are locally well described by two Fano forms. This, how¬ 
ever, is not surprising since by their very nature each of 
the separate resonances, one with a larger and one with 
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Table 2. Resonance energies, widths, and asymmetry param- Table 3. Fitting parameters for the cross sections shown in 
eters of the two Fano profiles shown in Fig. 0 Fig. [3] calculated using eq. 0- 


Ei r± E 2 r 2 



E 1 

A 

<li 

e 2 

r 2 

92 

a ) 

2.04970 

0.10792 

2.04990 

-0.09239 

a) 

2.0463 

0.0008 

-6.4770 

2.5416 

0.0995 

299.0411 

b) 

2.49559 

1.04493 

2.50780 

-0.91573 

b) 

2.0424 

0.0029 

-3.1026 

2.2997 

0.0969 

36.6894 

c) 

3.47188 

3.20257 

3.62857 

-2.28382 

c) 

2.0390 

0.0062 

-2.1014 

2.2028 

0.0934 

12.1206 

d) 

2.03912 

0.01471 

2.20748 

0.18673 

d) 

2.0244 

0.0438 

-1.8446 

2.0785 

0.0440 

1.8728 

e) 

2.04252 

0.00680 

2.30187 

0.19218 


f) 2.04626 0.00196 2.54172 0.19923 




energy E 


energy E 




energy E 


energy E 


Fig. 2. Fano fits of the cross sections shown in Fig[l] with two 
sets of Fano parameters. 


a smaller width, should be describable by Fano profiles, 
with different parameters. It is remarkable to see that the 
peaks retain Fano profiles even closer to and in particu¬ 
lar at the exceptional point, where the interference terms 
become increasingly important. 


3.3 Fit with energy dependent Fano parameters 


The idea of energy dependent Fano parameters has been 
discussed also in the context of q reversals in Rydberg 
series and multiphotoionisation 29 30 . Here we follow 
the formulation put forward by Magunov et al. [25] in 
their study of overlapping resonances. Starting from the 
S matrix with two isolated resonances Ek = Ek — ir k /2 
(k = 1,2), and a smooth reaction background described 
by a phase shift 26, they show that the cross section can 
be brought into the Beutler-Fano form 


o{E) = 4 sin 2 r/(E) 


(ei +q(E)y 

£? + l 


( 3 ) 


with 
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Fig. 3. Fano fits of cross sections at and in the vicinity of 
exceptional points using eq. 0. The blue curves are the cal¬ 
culated cross sections, the red curves denote the fits. 

The left-hand column shows the results for wi = 2.00 and 
LO 2 = 2.10,3.00,5.00 (from top to bottom) at the respective 
exceptional point (black dots), while in the right-hand column 
uj\ = 2.00, ui 2 = 2.10, and the distance from the exceptional 
point is successively increased, Af = 0.3,0.5,1 (from top to 
bottom). 


Apart from the scale factor for the height of the cross 
section, in the extended ansatz (0) the fitting parameters 
are E\, T), E 2 , A and the phase 6. 

In all our fits the numerical value of 6 turned out to be 
close to 7 r up to the first five digits. Since 6 enters into the 
S matrix by the factor exp(i2<5), this implies that in our 
model there is no background phase, and we could have 
chosen <5 = 0 from the outset. In this case 0 simplifies to 


a(E) = 


1 (ei+e 2 ) 2 

( e 2 + 1 ) (e? + l) 


( 4 ) 


r](E) = 6 - arccot(e 2 ), e k = 2(E - E k )/r k , 
and the energy dependent Fano parameter 
q(E) = — cot r)(E). 


The comparison with Q then shows that e 2 formally as¬ 
sumes the role of an energy dependent asymmetry param¬ 
eter q(E) for resonance 1, and vice versa. 

Examples for fits of the cross sections with energy de¬ 
pendent asymmetry parameters at and in the vicinity of 
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exceptional points are shown in Fig. [3] The agreement 
between the calculated cross section and the energy de¬ 
pendent Fano fit is found to be perfect over the complete 
energy range. 

The fitting parameters required for this excellent agree¬ 
ment are listed in Tab. 0 For sufficient distance from the 
EP (rows d-f), where the interference term plays a mi¬ 
nor role, this is to be expected as we encounter essentially 
two independent and well separated resonances associated 
with poles the positions of which are given in the table. 
In contrast, at the EP (rows a-c) the two peaks cannot 
be identified as two separated resonances; in fact, here 
the interference term is crucial to produce the two peaks 
(see 21 ). It therefore comes as no surprise that the best 
fits produce results - negative widths - that can no longer 
be interpreted as physical resonances. This clearly indi¬ 
cates that the treatment [25] by energy dependent Fano 
parameters fails at an EP. In fact, the double pole invoked 
by the EP is different in character from the double pole as 
considered in 125:: the two-dimensional coefficient matrix 
of the double pole at the EP does not have full rank but 
rank unity (see 20 ) meaning that the entries are corre¬ 


lated in a particular way. In this context, we also note that 
our scattering matrices are not unitary as the underlying 
Hamiltonian is not hermitian. 


4 Summary 

Starting from the fact that in scattering systems excep¬ 
tional points give rise to a second order pole in the Green’s 
function, we have investigated the effect on the shape of 
the scattering cross section as one approaches the excep¬ 
tional point. In this paper we have demonstrated that the 
asymmetric line profiles in the neighbourhood of the ex¬ 
ceptional point are to be interpreted as Fano resonances. 
Moreover, by allowing for energy dependent Fano para¬ 
meters the cross sections can also globally be described as 
actual Fano profiles. 

We feel that by this analysis we have given a new and 
deeper understanding of Fano resonances. On the one side 
there is the mathematical mechanism of two coalescing 
eigenvalues - an EP - on the other the physical origin of 
two interacting near resonances where the one can be a 
single particle resonance and the other a bound state in 
the continuum that, owing to the interaction, has acquired 
itself a width. This is the typical Fano-Feshbach scenario. 
Whether or not all resonances that can be fitted by the 
traditional Fano-Feshbach procedure have their origin in 
an EP cannot be shown. In fact, the particular type and 
relative strength of the interaction must be taken into ac¬ 
count. Yet we believe that the mechanism proposed here 
should be representative quite generally. 

It should be stressed that, in spite of the simplicity of 
the underlying classical problem, the generic behaviour of 
line profiles studied in this paper near exceptional points 
applies to any physical two-channel scattering or trans¬ 
mission problem. Therefore we did establish the direct link 
between the appearance of asymmetric Fano profiles and 
the occurrence of exceptional points in parameter space. 
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